Climbing the mountain: experimental design for the efficient optimization of stem cell bioprocessing

“To consult the statistician after an experiment is finished is often merely to ask him to conduct a post mortem examination. He can perhaps say what the experiment died of.” – R.A. Fisher While this idea is relevant across research scales, its importance becomes critical when dealing with the inherently large, complex and expensive process of preparing material for cell-based therapies (CBTs). Effective and economically viable CBTs will depend on the establishment of optimized protocols for the production of the necessary cell types. Our ability to do this will depend in turn on the capacity to efficiently search through a multi-dimensional problem space of possible protocols in a timely and cost-effective manner. In this review we discuss approaches to, and illustrate examples of the application of statistical design of experiments to stem cell bioprocess optimization.     

Molecular epidemiology and evolution of human respiratory syncytial virus and human metapneumovirus

Human respiratory syncytial virus (HRSV) and human metapneumovirus (HMPV) are ubiquitous respiratory pathogens of the Pneumovirinae subfamily of the Paramyxoviridae. Two major surface antigens are expressed by both viruses; the highly conserved fusion (F) protein, and the extremely diverse attachment (G) glycoprotein. Both viruses comprise two genetic groups, A and B. Circulation frequencies of the two genetic groups fluctuate for both viruses, giving rise to frequently observed switching of the predominantly circulating group. Nucleotide sequence data for the F and G gene regions of HRSV and HMPV variants from the UK, The Netherlands, Bangkok and data available from Genbank were used to identify clades of both viruses. Several contemporary circulating clades of HRSV and HMPV were identified by phylogenetic reconstructions. The molecular epidemiology and evolutionary dynamics of clades were modelled in parallel. Times of origin were determined and positively selected sites were identified. Sustained circulation of contemporary clades of both viruses for decades and their global dissemination demonstrated that switching of the predominant genetic group did not arise through the emergence of novel lineages each respiratory season, but through the fluctuating circulation frequencies of pre-existing lineages which undergo proliferative and eclipse phases. An abundance of sites were identified as positively selected within the G protein but not the F protein of both viruses. For HRSV, these were discordant with previously identified residues under selection, suggesting the virus can evade immune responses by generating diversity at multiple sites within linear epitopes. For both viruses, different sites were identified as positively selected between genetic groups.     

Evidence of Spread of the Emerging Infectious Disease, Finch Trichomonosis, by Migrating birds

Finch trichomonosis emerged in Great Britain in 2005 and led to epidemic mortality and a significant population decline of greenfinches, Carduelis chloris and chaffinches, Fringilla coelebs, in the central and western counties of England and Wales in the autumn of 2006. In this article, we show continued epidemic spread of the disease with a pronounced shift in geographical distribution towards eastern England in 2007. This was followed by international spread to southern Fennoscandia where cases were confirmed at multiple sites in the summer of 2008. Sequence data of the ITS1/5.8S/ITS2 ribosomal region and part of the small subunit (SSU) rRNA gene showed no variation between the British and Fennoscandian parasite strains of Trichomonas gallinae. Epidemiological and historical ring return data support bird migration as a plausible mechanism for the observed pattern of disease spread, and suggest the chaffinch as the most likely primary vector. This finding is novel since, although intuitive, confirmed disease spread by migratory birds is very rare and, when it has been recognised, this has generally been for diseases caused by viral pathogens. We believe this to be the first documented case of the spread of a protozoal emerging infectious disease by migrating birds.     

Cloning and heterologous expression of ectoine biosynthesis genes from Bacillus halodurans in Escherichia coli

The genes involved in the biosynthetic pathway of ectoine (2-methyl-1,4,5,6-tetrahydropyrimidine-4-carboxylic acid) from Bacillus halodurans were cloned as an operon and expressed in E. coli. Analysis of the deduced ectoine biosynthesis cluster amino acid sequence revealed that the ectoine operon contain 2,389 bp, encoded by three genes; ectA, ectB and ectC that encode proteins of 189, 427 and 129 amino acids with deduced molecular masses of 21,048, 47,120 and 14,797 Da respectively. Extracts of induced cells showed two bands at 41 kDa and 17 kDa, possibly corresponding to the products of the later two genes. However the expression of ectA gene could not be ascertained by SDS-PAGE. The activity of the ectA protein was confirmed by an acylation assay. The transgenic E. coli accumulated upto 4.6 mg ectoine/l culture. This is the first report of an engineered E. coli strain carrying the ectoine genes of the alkaliphilic bacterium, B. halodurans.     

Characterization and phylogenetic analysis of ectoine biosynthesis genes from <Emphasis Type="Italic">Bacillus halodurans</Emphasis>

Ectoine, a cyclic tetrahydropyrimidine (2-methyl-1,4,5,6-tetrahydropyrimidine-4-carboxylic acid), is a natural compound, which serves as a protective substance in many bacterial cells. In this study, the putative ectABC gene cluster from Bacillus halodurans was heterologously expressed in E. coli and the production of ectoine was confirmed by HPLC analysis. The activity of the enzymes coded by the ectA, B and C genes were found to be higher in induced transgenic cells compared to the uninduced cells. Phylogenetic analysis revealed sequence identities ranging from 36–73% for ectA gene, 55–81% for ectB gene and 55–80% for ectC gene indicating that the enzymes are evolutionarily well conserved.     

Genomic expression patterns distinguish long-term from short-term glioblastoma survivors: a preliminary feasibility study

We used microarray analysis to investigate associations between genotypic expression profiles and survival phenotypes in patients with primary glioblastoma (GBM). Tumor samples from 7 long-term glioblastoma survivors (>24 months) and 13 short-term survivors (<9 months) were analyzed to detect differential patterns of gene expression between these groups and to identify genotypic subclasses of glioblastomas that correlate with survival phenotypes. Five unsupervised and three supervised clustering algorithms consistently and accurately grouped the tumors into genotypic subgroups corresponding to the two clinical survival phenotypes. Three unique prospective mathematical classification algorithms were subsequently trained to use expression data to stratify unknown glioblastomas between survival groups and performed this task with 100% accuracy in validation studies. A set of 1478 genes with significant differential expression (p<0.01) between long-term and short-term survivors was identified, and additional mathematical filtering was used to isolate a 43-gene "fingerprint" that distinguished survival phenotypes. Differential regulation of a subset of these genes was confirmed using RT-PCR. Gene ontology analysis of the fingerprint demonstrated pathophysiologic functions for the gene products that are consistent with current models of tumor biology, suggesting that differential expression of these genes may contribute etiologically to the observed differences in survival. These results demonstrate that unique expression profiles characterize genotypic subsets of primary GBMs associated with differential survival phenotypes, and these profiles can be used in a prospective fashion to assign unknown tumors to survival groups. Future efforts will focus on building more robust classifiers and identifying additional subclasses of gliomas with phenotypic significance.     

Sparrowhawk Accipiter nisus presence and winter bird abundance

The decline and subsequent re-colonisation of the Sparrowhawk Accipiter nisus in the UK has allowed an examination of its potential impacts on passerine abundance through analysis of long-term garden bird monitoring data. Weekly counts in gardens in winter have been collected annually since the early 1970s, a period when Sparrowhawks were initially rare but recolonised to a point in the 1990s when they were recorded in the majority of gardens. Change in bird count from year-to-year for 10 common prey species from up to 483 survey sites were analysed in relation to Sparrowhawk re-colonisation, the number of feeding units per site (e.g. bird tables, hanging bird feeders) and, for a subset of the data, minimum temperature. Year-to-year change increased with Sparrowhawk re-colonisation in Great Tit Parus major, Blue Tit Cyanistes caeruleus and Collared Dove Stretopelia decaocto, but decreased in Starling Sturnus vulgaris, Chaffinch Fringilla coelebs and Greenfinch Carduelis chloris. Positive associations may suggest a behavioural response to predators, or an attraction of predators to increased prey caused by confounding effects of improvement in site quality. The significant negative effects detected suggest that there is a temporal matching between increased predation pressure and decreases in bird abundance for certain species. However, there was no strong evidence that Sparrowhawk re-colonisation was linked to long-term declines in population size.